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Abstract 

The distinction between the true total area and the projected area is elucidated with 
soluble models which represent the membrane as a self-avoiding string on a plane. Con- 
straining the total area to a predetermined value changes the averages very significantly. 
The latter are calculated exactly from the generating functions of self-avoiding walks and 
are shown as functions of activities q and r related to temperature T = ±1/ log(q) and 
lateral force / = — log(r). The constraint makes the partition functions and averages 
valid for all q, r > and reduces the ratio of A tot to the projected area L. High temper- 
ature divergences are supressed. Possible applications to simulated bilayers/membranes 

^ are discussed. 
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I. Introduction 

Interfaces and membranes as well all undergo shape fluctuations due to thermal 
motion. A two dimensional sheet such as a bilayer or a planar interface changes its 
shape by excursions into the third dimension; thus arises the distinction between the 
"true" total area At ru and the "projected" area. The latter term has been coined to 
emphasize that the "true" (proper, intrinsic) area of a planar interface when projected 
onto the edge of the volume, produces A, the projected, or nominal, area. The definition 
of the interfacial tension, via the const ant- volume increment 5F of the free energy F 

5F = ~f5A (1.1) 

refers to the projected area A . Another expression that both planar interfaces and 
membranes share, is the limiting form of the structure factor, dominated by capillary 
waves. 

However, there are important distinctions, essentially originating from the fact that 
an interface is an open system, also open with respect to particle exchange. A liquid 
surface fluctuates and varies its surface area by diffusion of molecules in and out of 
either bulk phase. A surface of a perfect crystal is an extreme example, as it varies 
exclusively by particle exchange with vapor phase. 

In membranes or bilayers such particle migration does not occur or is very rare. 
Commonly membranes and bilayers are formed by self-assembling surfactant (amphi- 
philic) molecules, embedded in a liquid solvent. These are practically insoluble in the 
solvent forming a variety of micelles if not in a bilayer. Consequently the shape fluctu- 
ations of a membrane or bilayer take place under the constraint of a constant particle 
number N s forming the surfactant sheet. 

This constraint translates to a constraint of constant membrane area A tru . 

The common formulation of the theory of capillary waves starts with an flat inter- 
face (with its true area A tru originally equal to the projected area) which increases its 
true area by thermal fluctuations. Thus the projected area is a given quantity whereas 
Atru varies in time about its thermal average (...) value. 

With a membrane one must take the opposite view: the true (intrinsic, proper) 
area of a membrane A tru is a given quantity whereas the projected area is a fluctuating 
quantity with a thermal average (...). 

The effects of the constraint of constant particle number in the membrane are 



discussed below. Introducing these constraints into the theory has proved difficult; 
practically the classical reference 1 ' 2 is entirely up to date 3 even today. A mention 
of some difficulties is present in a theorerical paper 4 . A most useful source of in- 
formation have been computer simulations. The areas of bilayers have also been ex- 
tracted from simulation data 5_7 . Overall, the simulations of bilayers have been very 
numerous 5-26 , but not all determined the bilayer isotherms i.e. lateral tension vs. area 
dependences 5 ' 6 ' 10 ' 12 ' 13 ' 16 ' 21 " 24 . 

Although the definition (1.1) of 7 may be the same, in fact the behaviour of 7 of 
bilayers has been so different from 7 of interfaces/surfaces - that a specific terminology 
of "lateral tension" was introduced. The "exotic" properties of the lateral tension have 
been recently enumerated 26 . 

In this paper we find useful insights from exact solutions of simple models. The 
membrane is modelled as a onedimensional string embedded in a square lattice. Fig.l 
shows such a string and the two areas associated with it. As calculations show, besides 
the the membrane/bilayer case, of chief interest, the interfacial situation can also be 
modelled with appropriate choice of parameters while removing the constraint. 

In Section II we apply the known generating functions 27 to derive some new variants 
of these and then to calculate the averages: the average projected area for given full 
area (the membrane) and the average total area for given projected area (an interface); 
we illustrate the spectacular differences with Figures. Finally we consider the case of 
a boxed membrane/bilayer, such as is modelled in molecular dynamics simulations, for 
which both areas, total and projected, are given and fixed. Section III is the summary 
and a discussion, illuminating i.a. the issue of the bending (rigidity) coefficient. 



II. The Random Walk and its Generating Functions. 

We model the membrane as a two dimensional string in an infinite square lattice. 
The string is generated by a partially directed self-avoiding random walk (PDSAW). 
With lattice constant ao = 1, the allowed steps are: up, down, or to the right, i.e. 
(0, +1), (0, —1), or (+1,0), respectively. The total number of steps, A, is made of L 
horizontal a;— steps, and of A v vertical steps. A v = n + + n_ is made of positive and 
negative steps n + > and n_ > 0. Thus always A = A v + L = n + + n_ + L. Such a 
string is shown in Fig.l. 

A string of length of A steps extends from (0,0) i.e. from x — 0, y — hi — 0, to 
(L,hi ast ) i.e. to x = L,y = fiL = hi ast . If the string is to represent a onedimensional 
membrane, the length A becomes the membrane area, and L the projected area. The 
usual definition of an area of a smooth surface translates into s 2 = Ax 2 + Ay 2 as the 
surface element produced by one step; we approximate the hypothenuse s by the sum 
s = Ax + \Ay\. Thus the area is always somewhat overestimated, but the difference 
with the usual definition was found to be quite small. Then the total area ("true") is 
A, Atm = A, and A = A v + L. In the membrane picture A once given is constant; it 
does not fluctuate. In the interface model, L is given and constant and A can take any 
value A > L. 

The generating function is constructed as 

all walks 

where r, q+, q- are weights assigned to these steps and counting parameters. Generally 
up (plus) steps are assigned the same weight as the down (minus) steps, i.e. q = q- = q + 
and q Av is substituted for q^_~ ■ q^ + . 

The parameter r > 1 will favor large L, i.e. extended, relatively flat configurations, 
whereas r < 1 will favor small L. The parameter q < 1 will favor less deviations from 
flatness, whereas q > 1 will favor configurations with many folds - which are suggestive of 
a crumpled membrane or floppy bilayer. Restriction of hi as t to the value hi as t = hi = 
may also be imposed in order to have periodic boundary condition in the x-direction. 

A still more general model would introduce the bending parameter w and a factor 
w T where T is the number of turns. Here we take w — 1, i.e. we ignore the number 
of turns. The issue of bending coefficient in the context of this model is discussed in 
Section III. 



Clearly (2.1) can be interpreted as a grand partition function with r, q as activities; 
thus 

q = exp[-/3e„] (2.2) 

where (3 = 1/kT and e v is the energy cost of a vertical step. Now we note that q < 1 
implies e v > i.e. the choice for modelling an interface which deviates from flatness 
at some cost. This cost is reflected in the positivity of 7 defined by (1.1). Otherwise, 
q > 1 and e v < owing to attractive intermolecular forces between the segments of the 
bilayer or membrane. 

The activity r may be interpreted as r = exp[— f3ch\ with some horizontal energy 
cost, but if we reintroduce the lattice constant ao, then we see that r L may be profitably 
interpreted as 

r L = exp[-/3/a L] (2.3) 

where now L is a nondimensional integer and / is a force in the x— direction. Then 
r > 1 favoring large L implies a pulling force (which is negative) whereas r < 1 implies 
a compressing force (which is positive); r = 1 is the locus of the tensionless states. 

We adapt to our purposes the general solution 27 for G 

G = (r + q-+q + - (r + 2)g_g+)/(l - r - q+ - q- + q-q+(l + r)) (2.4) 

(with w = 1 as mentioned above). Without distinction between positive and negative 
vertical steps, i.e. with q- = q + = q, 

G = (2q + r + qr)/(l — q — r — qr) (2.5) 

With the condition of at least one horizontal step, G is 

G\ — (t + qr)/(l — q — r — qr) (2-6) 

All three generating functions impose no restrictions on the final position (L, hi as t) 
attained after all iV steps; 

hi ast = n + - n_ + hi = n + - n_ . (2.7) 

For strings with periodic boundary condition, taken with Fourier analysis in view, the 
walk is from (0, 0) to (L, 0), h\ as t = . The appropriate generating functions with this 
constraint embedded, are derived below. 



The generating function G counts all horizontal L steps and all vertical A v steps 
which produce configurations of the string; the latter can be gathered together in the 
combinatorial factor g(L,A v ) 

G (q, r) = J2 Yl 9o(L, A v )q A ^r L N = A v + L (2.8) 

L A v 

and similarly G\. 

Canonical averages for an interface. 

The shape fluctuations of an interface take place at L fixed and constant. The 
statistical averages, e.g. (N) = (A v ) + L are calculated under this condition. From the 
generating functions, the Taylor expansion of the generating function 

G = J2 Z ( L ^)r L (2-9) 

L 

produces the partition function Z(L,q). The averages are found from 

{A v ) = q(d\ogZ(L,q)/dq) (2.10) 

and (N) = L + (A v ). The vertical part of N, A v , can be any nonnegative integer. In 
this way, for L > 2, from (2.5-2.6) 

Z(L,q) = «l + q)/(l-q)y L + 1 \ (2.11) 

The same Z obtains from G\. The average total area A^ follows from (2.10) 

(N) = L+(A v ) = L + 2qL/(l-q 2 ). (2.12) 

The areas are best represented by a normalized quantity C = L/N = L/(A V + L); 
always £ G [0, 1]. Rewriting (2.12) 

C = L/(L + (A v )) = / t (g) = (1 - q 2 )/(l + 2q- q 2 ). (2.13) 

This expression results from both G\ and G . Since L is constant, C is L/(N) = 
L/(L+(A V }). 

For the canonical averages at given q, r, i.e. without fixing L (or N) at a predeter- 
mined value, Gi gives fi(q,r), (2.13) again; However, from G by using 

(L) = r(dG /dr)/G (2.14) 
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one obtains 

C = oi(q, r) = (1 + q) 2 r/(2q + r + 2qr + q 2 r). (2.15) 

This does not differ much from (2.13) and is quoted only for completeness. The inde- 
pendent variables were L, q or q, r; the length N and A v resulted as averages, so that C 

was (L)/((L} + (A V ). 

Keeping the distinction between q + and q- and expanding (2.1) in powers of r, we 
find Z(L, q + , </_) = \1/ L where 

1-q 2 
1 - Q+ ~ Q- + Q 2 
and the expansion of Z is 



#=- : q ■ 2 . (2.16) 



Z(L,q + ,q-) =^2g(n + ,n-,L)q r l + q 1 l (2.17) 

The walk beginning at (0,0) and ending at (L, hi ast ) will have fa/ aS i = hi ast — h\ — 

n+ — n_. Let us call this quantity n. Then in the sum we insert a Kronecker Delta 

ensuring that n = n+ — n_; by using its Fourier representation, putting q + = q- = q, 

we obtain 

Z(n;L,q) = (l/ir) f dkcos(n* k)( 1 " q 2 ) m (2-18) 

Jo l — zq cos{K) -f- q 

where the power misL, L + l, L — 1 depending on the starting generating function. The 

result of integration depends on the inequality (1 + q 2 ) > \2q\ which is fulfilled for all 

q > except for one point q = 1. This point approached from below corresponds to the 

infinite temperature in the interface models with q = exp(— (3e v ), e v > 0; approached 

from above corresponds to the infinite temperature in membrane models for which 

e v < 0. For n = in (2.18) we obtain the partition function for an interface with 

periodic boundary condition; it reads 

Z(0;L,g) = P m _ 1 (u) u = (1 + q 2 )/(y/l - 2q 2 + q 4 ) (2.19) 

Here P m (u) is the Legendre polynomial, u > 1. The normalized average area C is in 
this case L/(L+ (A v )) and (A v ) is calculated after (2.10) as q(d log Z/dq); we find then 
that we must take q < 1. Finally 

C=(l + u-P L - 1 (u)/P L (u))- 1 . (2.20) 

As we shall see below, this quantity behaves quite unlike those £'s constrained by the 
total length N = L + A v . It agrees very closely with fi(q, r). 



Canonical averages for a membrane/bilayer. 

Alternatively, treating the string as a membrane, we impose a fixed length N. This 
constraint is introduced into G's. The new generating function contains a selection of 
terms from (2.5), such that L + A v = N = const., 

N 

G(N;q,r) = ^g{A v = N - L,L)q N - L r L N = A v + L = const. (2.21) 

L=l 

There are two cases possible: either all L compatible with given N are allowed, or both 

L and N are given prescribed values. In the latter case the expansion (2.21) reduces 

to one term, with (A v = N — L,L). From the expansion (2.8) of G , in powers of r 

and then q, we find g (A v , L); alleviating the notation by writing a = A v , we write the 

explicit result 

v^ fL\fa — m + L — l\ . n . 

*("•*> = E(J( L-i ) < 2 - 22 ' 

m=0 V 7 V 7 

For a > L the upper limit of the sum is L. With g known, we can construct the 
partition function Z as a polynomial in q and r, 

JV 

Z {N- q,r) = Y, 9o(N -L,N)x q^^ x r L (2.23) 

by explicit substitution of (2.22) into (2.23). Now (L) and (a) follow explicitely by 
(2.10) and (2.14). The same procedure works for G\. For these averages at q,r;N 
always C = (L)/((L) + (Av)). 

A direct and neat way is to introduce a new counting parameter t for every step 
whatever it is and expand 

G (t; q, r) = (2qt + rt + qrt 2 )/(l - qt - rt - qrt 2 ) (2.24) 

in powers of t, 

oo 

G (t; q,r) = J2 Z n(q, r)t N (N = A V + L). (2.25) 

t-i 

Zn is used for calculations of average (L) (which must be equal to A^ — (A v ). 

Finally, pinning both ends of the string in order to have a possiblility of p.b.c, is 
done in two alternative ways: either we expand (2.4) in powers oft, select the coefficient 



of t N which is Z^(q+, q-,r); substitute there g_ — > q/e and q + — > qe, expand in powers 
of q and e and select the coefficient of e° - it is the partition function Z^(q, r) for both 
ends pinned, useful for p.b.c. Alternatively, the equality of n + and n_ is enforced via a 
Kronecker Delta and its Fourier representation as an integral. Either way, its generating 
function is found with a rather unusual appearance 

G e (q,r) = r*VB;B= (1 - q 2 )/((l - q - r - qr)(l + q - r - qr)). (2.26) 

To calculate the averages at constant N, the same procedure is used as for G . Put 
q — ► qt,r — ► rt, extract the Taylor series coefficient of t N and use (2. 10), (2. 14). The 
resulting le(q,r) depends on N. 

These exact enumerations can be pushed very far to iV > 100 - but then simple 
asymptotic expressions take over. These are derived now. 

Both G\ and G are positive and finite for small q and r; their divergences are 
determined by zeros of the denominators. This condition, 1 — r — q — qr = 0, determines 
a line in the q, r plane cutting diagonally from (0, 1) to (1, 0). We write the denominator 
as a(t — t\){t — t 2 ) where t\^ are its roots and split G\ as 

G 1 =B 1 /(t-t 1 ) + B 2 /(t-t 2 ) (2.27) 

The root t\ is always positive and t\ < 1/r; the other root is always negative. G 
takes largest values when t is close to t\\ expanding in powers of t (as (1 — £/£i)( -1 ) = 
1 + (t/ti) + ■ ■ ■) we see that the root t\ > is the one that matters. From this new form 
of G\ we select the coefficient of t N and calculate (L), taking the asymptotic limit of 
large power of t, i.e. large N. Finally we obtain 

C = f 2 (q, r) = (1/2) + (l/2)(r - q)/^{r 2 + 6rq + q 2 ) N = const. -> oo. (2.28) 

This expression, like (2.22-2.24), is valid in the entire quarter-plane (q > 0, r > 0). £ 
from G attains the same asymptotic result in the limit iV — ► oo. 

An important point is that the partition functions for fixed N (total area, or length 
of the string) Z^(r,q) are positive and well-behaved for all q > 0, r > whereas the 
generating functions G(r, q), the partition functions for fixed L, Z(L, q) and the averages 
and fluctuations derived therefrom, are well behaved and physically acceptable only in 
certain regions of the quarter-plane. 
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Fig. 2 shows five functions C: (1) the canonical average at given q, r, A, denoted 
as ll(q,r); eq. (2. 25), (2. 10), (2. 14); (2) as (1) but under the periodic boundary condition 
(i.e. pinning both ends at hi ast = hi) - function le(q,r), eq.(2.26)ff; (3) the asymptotic 
limit of (1) for A — *■ oo, denoted as f2(q, r), eq.(2.27). (4) the canonical average at fixed 
L (and given q), any A, any hi as u function fi(q,r), eq.(2.13). (5) as (4) but under the 
periodic boundary condition hi as t = h±; eq.(2.20). 

In order to be physically acceptable, each £ must be C G [0, 1] and also the averages 
(L) and (A v ) must be positive. This selection of 5 functions £ are shown plotted against 
q, for a series of values or r. The regions q < 1 and q > 1 are smoothly joined by 
all averages calculated at constant A. The function f\ behaves quite differently: its 
physical branch is limited to the region < q < 1. f± does represent an interface; the 
other three - a membrane. Similarly under the periodic boundary condition C(q,L), 
(eq.2.20), for either L = 10 or L = 140 does not differ visibly (on the scale of the 
plot) from /i; it also represents an interface. All functions start from A v = at q = 
but the interface fluctuations give unbounded (A v ) at q — *■ 1~ whereas the membrane 
fluctuations are bounded. The three constant-iV functions ll(q,r;N),le(q,r;N) and 
f2(q, v) cross smoothly the point q = 1; all three are not too far away from each other. 
In particular the asymptotic J2 would be merged with the other two on the scale of the 
Figures if the plots were not drawn for a choice of small A. 

The interval < q < oo is split in two and in Fig. 3 amd 4 the plots are against 
inverse temperature, against b = —log(q) for q < 1; against b = log(q) for q > 1, 
respectively. In Fig. 3 the limiting values for 6^0 (the limit of high temperatures) 
depend on the parameter r, except for the two cases of fixed projected area L where 
the limit is C = 0, as L=const. and (A v ) is unbounded. The point b = + corresponds 
to q — *• 1~. f2(q = l,r) (eq.2.27) is continuous and /2(1, 1) = 1/2. All curves, fixed 
A^ or fixed L, merge towards C = 1 as b increases, i.e. kT — *■ 0. Obviously, (A v ) 
vanishes there. In Fig. 4 the function /i is shown and how it takes unphysical values 
either negative or larger than 1. The other L— fix average(eq.(2.20) is also unphysical in 
this region, unlike the three constant-A functions. At very low temperatures all three 
£'s tend to zero because A v dominates; q ^> 1 favors configurations with many folds for 
which A v 3> L. 

We have chosen to convert all averages to £'s for convenience but in this way some 
dramatic differences in behaviour are hidden from view. Thus C = implies unbounded 
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A v at finite L or vanishing (L) at finite A v ; C = 1 if A v = 0. In the limit of high 
temperatures, the interfacial averages, i.e. those taken at constant L, diverge, whereas 
the membrane averages i.e. those at constant N, produce finite values (depending on 

r). 

Next we change the role of variables and plot below in Fig. 5 the same functions 
against r at several constant values of q; the logarithic scale implies that we choose force 
/ = — log(r) as the independent variable (the ordinate). 

A boxed membrane. 

In simulations most often not only the number of molecules forming the bilay- 
er /membrane is kept constant, but also the membrane or bilayer is confined to the 
simulation box of dimensions L x x L y x L z . For our model this implies that not only 
N is given and constant but also L. Because in this particular model the sum A v + L 
coincides with N, at any given temperature the partition function and the free energy 
are 

Z = g(A v ,L)q A " (A V = N -L,N = const., L = const.) (2.29) 

[3F = -logg{A v ,L)-A v logq (2.30) 

The activity q = exp[e//cT] > 1. For comparison, a boxed interface in an identical 
simulation box, does not keep A v + L constant and therefore 

Zint = J2 9 ^ ^^ ( L = COnSt ^ ( 2 - 31 ) 

a 

[3F int = -logZ tnt (2.32) 

The lateral tension of the membrane results as 

OF 
T = {-qj)t,n = -{l/g){dg/dL) N (2.33) 

The derivative dg/dL of g(A v ,L) is taken at constant iV . Fig. 5 shows C for a series 
of g's with constant iV plotted against the negative of the force (— /) = log(r). The 
asymptotic f2(q,r) (for iV — * oo) is also included - it has very much the same shape, 
only shifted along the ordinate according to the value of q, just like all the finite— N 
curves. For comparison also T from finite differences (Ag / AL) n is shown demonstrat- 
ing the equivalence of ensembles. The continuity w.r.to q is also seen. The pulling 
force raises (L) as expected and conversely; the effect of q is to shift the entire curve 
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L(f) to the right or to the left, without much of a visible change of shape. This can 
be understood by looking at the maximum term in the canonical sum for which the 
condition is (dlogg/dL)^ = —f — log q. The parameter q shifts the curve f*(L) by a 
constant amount. 

III. Summary and Discussion. 

The model allows explicit calculations in terms of simple algebra. Its use made 
it possible to demonstrate the role of the membrane-constraint in the statistical me- 
chanics of membranes or/and bilayers. The constraint of constant particle number in 
the membrane - no particle loss or gain to or from the surrounding solvent - calls for 
a corresponding approach in the calculation of the membrane thermal fluctuations. In 
this model this translates into a constant length of the string. Such a constraint is most 
natural for modelling of a polymer and we expect that the calculations for a twodi- 
mensional membrane embedded in three dimensions, will be in effect calculations for a 
two-dimensional polymer. 

In the course of calculations we have found that relaxing the constraint of constant 
N and replacing it by a constraint of a fixed projection area L, leads to a very different 
behaviour of the averages and - a posteriori obviously - modelling the fluctuations of an 
interface. Thus in the same model we have the region q > 1 inaccessible to interfaces but 
typical for a membrane - and the region q < 1 typical for interfaces and also accessible 
for a membrane with a changed sign of the energy of interaction. 

The equivalence of ensembles operates fully in the limit of infinite systems but not 
for finite sizes. 

The force, pulling or compressing the string, has a direct analogy in three dimen- 
sions, as the lateral tension of the bilayer. Cutting a suitably chosen part of an S-shaped 
curve from Fig. 5, one can get a picture very similar to the "bilayer isotherm" (which 
is a plot of lateral tension against the projected area). Curiously, a non-symmetric 
derivative (dg/dL)^ i.e. (g(L + 1) — g(L ))/g(L ) assigned to the point L , produces 
such curves. 

Because the string is confined to the square lattice, it has an implicit resistance 
to bending which is also seen in the play of the parameter w. The latter has been 
introduced in the generating functions of the PDSAW 27 in the form of w T where T is 
the number of turns; see also 28-30 . Clearly w ^> 1 will favor the occurence of multifolded 
configurations with many turns - very low bending coefficient, low resistance to bending. 
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Conversely, w — > will favor "no turns" i.e. flat membranes which are obtained with 
large bending coefficients - large resistance to bending. Thus it follows that our choice 
of w = 1 does not mean "zero bending coefficient" - and the explanation lies in the 
implicit resistance to bending imposed by the lattice. The notions of " semiflexible" and 
"super-flexible" strings have been introduced 28-30 . 

In this context it becones necessary to mention the existence of a sizable amount of 
very successful work on exact enumeration - for models of polymers as strings confined 
to a lattice 28-34 . References to earlier work can be found there as well. However, the 
issues tackled in this paper were not touched in those References, where the emphasis 
was entirely on the important topics of phase transitions, singularities, and critical 
scalings. 

Although the usual picture of a membrane assumes an attractive interaction be- 
tween its segments - which leads to q > 1, it is possible to envisage a string held together 
by other (intramolecular) forces while the segments repel each other or a good solvent 
encourages extended configurations; this leads to q < 1 - still a membrane if N is held 
constant. 

Incidentally we also confirm the view 26 of the tensionless state / = 0, r = I as 
just another point on the r axis or line on the q, r plane without any special features. 
Contrariwise, q = 1 is a very special point as discussed in Section II and shown in the 
Figures. 
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Figure Captions. 

Caption to Fig.l 

The string produced by a random walk on the square lattice. The x coordinate 
measures the projected distance-area L, the y-coordinate is called "height" in the text. 
Here the walk starts at hi = and ends at hi as t — hi. In the Figure the number of 
"horizontal" steps is 14, n + = ri- = A v /2 =, and the walk is a PDSAW. 

Caption to Fig. 2 

Normalized £, (L)/N, (L)/(N) or L/(N), as function of q, r, plotted against q G 
[0, +oo] for several values of r. Thick solid line: function f\ eq.(2.13) - unrestricted 
canonical average at given (q, r) coincidentally equal to average at fixed L. Diamonds: 
fixed L = 140 and pinned at h\ = hi ast = ( for the p.b.c.);eq.(2.20). Membrane 
averages for fixed N: thin broken lines - asymptotic function f2(q,f) eq.(2.27); thick 
broken lines - function le(r,q), eq.(2.26) and ff. text; thin lines - function ll(q,r) - 
canonical average under the constraint of imposed value of N; here a small value of 
A=10 is chosen to emphasize the differences. All three £'s are moving smoothly with 
change of the parameter r; the Figure shows r = 0.1, 1., 9., 30. from lowest values to 
highest. See text. 

Caption to Fig. 3 

Normalized (L)/N or (L)/(N) or L/(N) plotted against inverse temperature b = 
— \ogq = e/kT for < q < 1. See Caption to Fig. 2. Thick solid line: function /i 
eq.(2.13), fixed L; Diamonds: fixed L = 140 and p.b.c. (hi as t = 0), eq.(2.20); the 
differences with /i cannot be seen on the scale of this graph. Thin broken lines - 
asymptotic function f2(q,r) eq.(2.27); thick broken lines - function le(r,q), eq.(2.26) 
and ff. text; thin lines - function ll(q,r) both for a small value of iV=10. The high- 
temperature limits depend on r, q, except for fixed L averages where the limit is zero. 
The low temperature limit is unity for all five £'s. 

Caption to Fig. 4 

Normalized (L)/N or (L)/(N) or L/{N) plotted against inverse temperature b = 
+ \ogq = e/kT for 1 < q < oo. see Caption to Fig. 2, 3. Thick solid line: function f\ 
eq.(2.13), fixed L ; all values unphysical. The other L at constant A^: with thin broken 
lines - asymptotic function f2(q,r) eq.(2.27); thick broken lines - function le(r,q) for 
pinned string, eq.(2.26) and ff. text; thin lines - function ll(q,r) for floating string - 
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both for a small value of N=10. The high-temperature limits vary with r, q. ?The low 
temperature limit is unity for all three averages under constant N. ?The low tempera- 
ture limit is zero, as q ^> 1 favors multifolded configurations with (A v ) dominating over 
(L). See text. 

Caption to Fig. 5 

Plot of £ at selected values of q against r. The logarithmic scale for r produces the 
negative of force -/ = log[r]. N = 40, q = 0.01, 0.1, 0.5, 0.9, 1, 100, 10 4 (thin lines in 
the sequence from left to right). Change of q shifts without change in shape (see text). 
One asymptotic f2(Q,r), eq.(2.27), is shown with a thick line. Diamonds and crosses: 
force from eq.(2.33) for boxed interface with N = 40 and N = 60; triangles for N = 40 
with a different interpolation of (1 / g)(dg j 'dL) n . Also see text. 
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